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(UNCUSSIFIED)  ABSTRACT 

This  report  explains  how  one  may  use  the  AFAMF  code  to  compute 
approximately  the  flow  arising  from  the  detonation  of  a cloud  of  fuel 
air  mixture.  A scaling  method  is  given  to  effect  a comparison  of 
computed  results  with  experimentally  measured  values.  ^ 
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LIST  OF  SYMBOLS 


radial  distance 

r.  dimensionless  radial  distance 
J 

radius  of  volume  equivalent  spherical  cloud 

Pj  detonation  pressure 

pressure  at  the  surface  of  the  cloud  (or  piston) 

P^  ambient  atmosphere  pressure 

Pj^  shock  overpressure 

2N  number  of  mesh  points  in  initial  wave 

2.  dimensionless  variable 
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(UNCUSSIFIED) 


1. 


INTRODUCTION 


In  an  earlier  report  [I]  a model  has  been  developed  to  approximate 
the  flow  which  arises  from  the  detonation  of  a spherical  cloud  of  a 
fuel-air  mixture.  This  model  uses  a G.  I.  Taylor  [2]  self  similar  wave  to 
simulate  the  flow  during  the  detonation  process.  After  the  detonation 
ends  the  flow  is  computed  by  finite  difference  methods.  The  surface  of 
the  cloud  is  held  stationary  after  the  end  of  the  detonation  process  to 
avoid  the  computational  difficulties  of  a contact  surface.  This  imposed 
constraint  has  the  effect  of  starting  a rarefaction  wave  propagating  into 
the  flow  field  in  our  computation  earlier  than  it  happens  in  the  actual 
physical  phenomena.  In  the  earlier  report  [1]  it  is  shown  how  one  computes 
various  quantities  which  are  important  in  the  assessment  of  the  damage 
resulting  from  such  clouds.  Experimental  results  exist  for  a variety  of 
clouds  of  various  shapes  and  of  various  fuels.  In  this  report  it  is 
intended  to  show  how  one  can  scale  these  results  in  order  to  effect  a 
comparison  with  the  computations  given  in  [1]. 

(UNCLASSIFIED)  II.  SCALI.NG  OF  RESULTS  AND 
COMPARISON  WITH  EXPERIMENT 

Given  the  chemical  composition  of  the  cloud  of  fuel  and  the  mole 
percentage  of  fuel  to  air  in  the  mixture  one  can  compute  the  detonation 
pressure  P^  of  this  mixture.  To  compute  P^  there  are  a variety  of  codes 
which  can  do  this  such  as  TIGER,  TAMER,  and  FAE.  Through  an  interpolation 
one  determines  a value  of  a which  gives  a wave  with  the  pressure  at  the 
surface  of  the  cloud  P^  equal  to  P^.  Table  I gives  corresponding  values 
of  a and  Pg/Pg  which  can  be  usefully  employed  for  such  a purpose. 
Essentially  the  parameter  a determines  the  piston's  velocity  which 
produces  the  desired  wave.  Since  most  of  the  available  experimental 
results  are  for  clouds  which  are  not  spherical  one  has  to  measure  the 
volume  of  the  cloud  and  then  compute  the  radius  R^  of  the  sphere  which 
has  an  equal  voliome.  From  the  computations  one  obtains  the  shock 


*References  are  listed  on  page  12. 
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over  pressure  as  a function  of  the  position  r^ 
distance,  where 


in  dimensionless  radial 


Tj  - i ♦ (j-2)Ar  , 

4r  » (Exp(:j^)  - 1)/2N  . 
To  convert  to  dimensional  distances  let 


(1) 

(2) 

(3) 


A plot  of  Vs.  r^  is  made  to  compare  with  those  measured  experimentally 
after  converting  P.  to  Psi.  For  a 4’,  mixture  of  CH,CH,CH,NO,  normal-propyl 
nitrate  (NPN)  and  96%  air  it  was  found  that  the  detonation  pressure 
Pj  = 18  atmospheres.  The  interpolated  value  of  a « 1.77616  gives 

= 17.9865.  The  flow  computations  of  this  case  was  performed  using 
the  AFAMF  code  developed  and  described  in  [1].  A plot  of  the  computed 
shock  overpressure  as  a function  of  distance  was  made  in  Figure  1 for  a 
spherical  cloud  of  radius  » 11.6  ft.  of  the  above  described  composition. 
Experimental  data  were  later  on  obtained  from  TBL  and  plotted  in  Figure  1 
to  effect  a comparison.  The  experimental  values  plotted  in  Figure  1 are 
mean  values  in  the  sense  that  they  are  taken  from  a curvilinear  fit  of 
mean  experimental  values.  Figure  2 shows  computed  values  of  the  static 
and  dynamic  impulse  for  such  a cloud  plotted  as  functions  of  distance  from 
the  center  of  the  cloud.  The  comparison  of  experimental  and  computed 
overpressures  indicates  that  this  code  can  be  profitably  used  to  predict 
the  flow  from  such  clouds.  From  Figure  1 it  appears  that  the  approximation 
is  a good  one.  To  compute  the  reflected  shock  pressures  one  can  refer 
to  the  theory  developed  in  [3],  or  for  normal  reflections  use  the 
results  of  Shear  and  McCane  [4]. 
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TABLE  I (U) 

a 

P /P 

S 0 

a 

a 

%/'’o 

0 

1.0 

.8 

2.1389 

1.6 

10.3751 

.1 

1.0189 

.9 

2.5108 

1.7 

13.8711 

-> 

1.0704 

1.0 

2.9725 

1.8 

19.7021 

.3 

1.1517 

1.1 

3.5482 

1.9 

31.2932 

.4 

1.2639 

1.2 

4.2734 

2.0 

65.2023 

.5 

1.4113 

1.3 

5.2012 

2.1 

1676.1990 

.6 

1.6004 

1.4 

6.4157 

.7 

1.8397 

1.5 

8.0566 
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